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Abstract. Clumping in stellar winds of hot stars is a possible consequence of radiati- 
ve-acoustic instability appearing in solutions of radiative-hydrodynamical equations. 
However, clumping is usually included to stellar atmosphere modeling and radiative 
transfer calculations in a highly approximate way via a global free parameter called 
the clumping factor Using different values of clumping factors many researchers suc- 
ceeded to fit the observed spectra better and to correct empirical mass loss rates. This 
i-^ ■ usually leads to a conclusion that the stellar wind is clumped. To understand how 

clumping may influence theoretical predictions of mass-loss rates, diff'erent clumping 
properties have to be taken into account. If clumping appears already below the critical 
point, the mass-loss rates is changed. 



1. Introduction 

Many inhomogeneous structures in different astronomical objects have been directly 
observed, in supernova remnants, pla netary nebulae, or ju st recently structures in the 
envelope of the red supergiant a Ori (|Kervella et al.ll201lh . Consequently, there is no 
reason to assume that the winds of hot stars are perfectly homogeneous. 

However, there are no such direct observations of wind structures available. The 



presence of clumping w as inferred from lin e profile variability (e.g. Eversberg et al. 
19981 : iLepine &Mofeill999l : iMarkova et aDl2005l : |Hamann et al "2008'). Indirect ev- 



idence of clumping comes from X-ray observations dOskinova et al. 2004,) . Another 
evidence of clumping follows from detection of X-ray flaring in V ela X-1, caused prob- 



irmg in v et 

ably by wind clumps falling onto a neutron star (iFilrst et al.ll2010l) . 

The generally accepted origin of clumping in hot star winds is the radiative-acoust- 
ic instability, supported by several numerical simulations (see the review by Sundqvist, 
these proceedings). An alternative explanation using adiabatic fluctuations (IChiuehl 
Il997l) did not attract attention of research ers. However, recent interesting idea of the 
influence of subphotospheric convection dCantiello et al.l 12009) offers a possibility of 
existence of clumps in the part of the w ind below the critical point, where the numerical 
wind simulations do not predict it (e.g. lFeldmeieilll995l : lRunacres & Owockii.200 2). It 



seems that the most probable clumping scenario consists of multiple clump formation 
mechanisms. Although there is a lot of unknown about clump formation, distribution, 
and properties, to find some kind of description of clumping is very important. 
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2. Description of clumping 
2.1. Void interclump medium 

Due to missing direct observation of wind clumps the detailed treatment of clumping is 
usually restricted to using free parameters. This approach was previously used for the 
description of condensations in nebulae by lOsterbrock & Flatbed (11959 ^ for the case 



of the Orion nebula, where they assumed void intercondensation medium and used 
the size of condensations as the free parameter, since the observational data available to 
them did not allow to resolve details. Contemporary advanced observational techniques 
made the usage of the free parameter unnecessary. 

The assumption of the void interclump medium is a common one for the descrip- 
tion of a clumped wind. In other words it means that all matter is concentrated to 
clumps in vacuum. Such medium can be described by means of one free parameter, 
which may be either the (volume) filling factor 

y _ ^clumps 
^wind 

a fractiona l volume which contains material at higher density (referred to as the filling 
fraction by lOwocki & Cohenll2006l) . or the clumping (correction) factor 



(Pwind) / 

Just one of these two adjustable parameters is sufficient. However, there is no reason to 
assume that the parameter defined by ([Hi or © is depth indepen dent. Its depth depen - 
dence follow s already from the pioneering hydrosimulations of lOwocki et al.l ([IQSS'). 



Hillier & Miller, (,1999i) introduced an expression for the depth dependence of the filling 



factor O 



/W = /oo + (l-/oo)exp|-^|, (3) 



where fco corresponds to ([T|l and t^ci is the location in the wind where clumping becomes 
important. Many radiative transfer calculations have been done using this expression. 
Recently, Puis et al. (,2006.) tried to determine the depth dependence of the filhng factor 
from observations. 



2.2. Dense interclump medium 

The assumption of the void interclump medium simplifies the description of clump- 
ing. The non-void interclump medium was assumed already by lAbbott et al.l (Il98lh . 
however, it was abando ned by later s t udies. Recently the idea of a dense interclump 
medium was revoked by Zsargo et al. ( 20081). The ef fect of dense interclump medium 
was studied in detail by Sundqvist et al. ( 2010ll20Tll) and Surlan et al. (these proceed- 
ings). Inclusion of the non-void interclump medium requires additional free parameter 
d, which relates the interclump density to the clump density or to the density of the 
smooth wind, respectively. 



Validity of clumping approximations for mass-loss rates determination 



3 



2.3. Clump properties 

In most radiative transfer calculations with clumping it is assumed that clumps are 
optically thin, which means that clumps are smaller than the mean free path of photons. 
This assumption is sometimes referred to as the "microclumping". However, it is more 
natural to assume that clumps may be optically thick in some frequencies. These clumps 
are larger than the mean free path of photons, which may happen both in co ntinua 
and li nes. This situa tion is being referered to as the "macroclumping" (Oskinova et alj 
20071) or "porosity" (iQwocki & Cohenll2006l) . 

Since no 3D hydrodynamical simulations of the wind are available, there is no 
hint what the shape of clumps may be. In the parametric treatment using the clumping 
factor nothing is explicitly assumed about the clump shape and sizes. In more detailed 
calculations different shapes were assumed, like spheres ( Surlan et al., these proceed- 
ings), cubes (.Muijres et al..,2011) . shell fragments (Oskinova et al.ll2004l) . On the other 
hand, it is not clear how important is the detailed clump shape, probably more important 
factor is the distribution of clumps due to the stochastic nature of clumping. 

The common assumption of most calculations is a smooth velocity field both in- 
side and outside clumps. However, as follows from hydrodynamical simulat ions, as 



suming nonmonotonic field inside clumps is probably closer to reality (see IQwocki 
20081) . This generalization was studied in more detail by ISundqvist et"aLl (l2010Ll201lf 



and these proceedings) and Surlan et al. (these proceedings). Note, however, that inho- 
mogeneous velocity field can affect only spectral lines. 



3. Influence of clumping on empirical mass-loss rates 

Mass-loss rates for particular stars are usually determined from comparison of model 
emergent radiation with observations. The emergent radiation is usually calculated 
assuming (i.e. for given) velocity v{r) and density p(r) structure (and, consequently, 
the mass-loss rate M and the terminal velocity Va,)- The velocity structure is usually 
assumed to obey the so-called yS-velocity law v = Vcoi^ - Rt/r'f, where is the stellar 
radius and yS is a free parame ter. Note that th is type of dependence was already derived 
by Milne ( 1926) and IChandrasek hai] (Il934l) with = 0.5. Today's estimates of this 



value are a bit higher. 

The mass-loss rates are usually determined using different parts of the stellar 
spectrum, like radio flux, infrared flux, the Ha line, UV resonance lines, and using 
synthetic s pectra from mod el atmospheres. Various determination methods were re- 
viewed by Puis et al. ('2008'). Different diagnostics result in different mass-loss rates 
(e.g. Bouret et al. 2003 ; Fullerton et al . 2006). Clumping in stellar winds is the sug- 
gest ed and promis ing way out from this problem. 

jHiUier ( 199ll) tested clumping by artificial periodic variation of p(r). Later lHillier & MiUei 



(Il999h introduced a depth variable filling factor which became a common method 
for inclusion of clumping in line formation calculations. Since many studies found that 
mass-loss rates with clumping taken into account in this way fit observations better, 
then the conclusion about the presence of clumping was drawn. 

3.1. Unsaturated resonance lines 



Saturated resonance ultraviolet lines fit theoretical spectra usually well, but their sen- 
sitivity to changes of wind parameters is relatively low. On the other hand, if lines are 
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not saturated, they become a sensitive diagnostics tool. This happens for resonance 
lines of less abundant ions like S iv, S v, Si iv, or P y. The case of the Si iv resonance 
hne (1393.76 A, 1402.77 A) was recently studied bv lPrinja & Massal dlOlOh and from 
similar line strengths of the resonance doublet they concluded that the wind is clumped, 
since in the opposite case the line ratio should be 2. 

The case of the resonance lines of P v (1 1 17.98 A, 1 128.01 A) is more famous. The 
po ssible importance of th i s doublet for mass-loss rates determination was po inted out 
by lCrowther et all (|2002| ). iMassa et all (l2003l ) and later iFullerton et all (12006 ) found a 
discrepancy in mass-loss rates determination and concluded that either the mass-loss 
rates determined from P v resonance lines are wrong (which may be corrected by in- 
clusion of wind clumping) or the a bundance of P v is lower . 

To test the abundance of P v, iKrticka & Kubat (120091) studied the NLTE ioniza- 
tion balance of phosphorus in the wind. Using the code of iKrticka & Kubali (120041) 
they calculated the phosphorus ionization balance both with and without X-rays and 
found that the changes were insignificant. Consequently, the abundance of P v is not 
lowered by presence of X-r ays, which supports the clumping hypothesis. Howerver, 
IWaldron & Cassinellil (l2010l ) suggested that XUV radiation may be important for phos- 
phorus ionization balance, detailed NLTE calculation is still missing, however. 



4. Influence of clumping on predicted mass-loss rates 

While it is basically known how clumping (treated in the parametric way using the 
volume filling factor) influences the empirical determination of mass-loss rate for given 
velocity and density structure, for the case of mass-loss rate predictions the situation is 
much less clear. 

Solution of hydrodynamic equations for given basic stellar parameters (e.g. the 
effective temperature Teff, luminosity L*, and radius /?*) and radiation at the lower 
boundary gives the density p{r) and velocity v(r) structure of the wind. The values 
of the mass loss rate M and terminal velocity Voo follow directly from the determined 
structure, they are usually called predicted. The principal question is, how clumping 
influences predicted value s of M and Vop . 

This was studied by 'Krticka et al.' ('2008a'.'b^ using the stationary hydrodynamic 
code for NLTE stellar winds ( Krticka & Kubat 2004 ). which consistently calculates 
the radiative force without using the parameterization by the CAK parameters k, a, 
and 6). Clumping was treated in an approximate way (pdump = Cdp) - cf. Eq.jJJl, void 
interclump medium was assumed. This assumption modified the opacity and emissivity 
in a different way for lines and continua, consequently the radiative force was modified. 
Both optically thin and optically thick clumps were tested. 

They found that if clumping starts below the critical point, then the mass-loss 
rate increases. If clumping starts above the critical point, then the mass-loss rate does 
not change, but the terminal velocity increases. Clumps larger than the Sobolev length 
result in the decrease of the mass-loss rate. Clumping also influences the ionization bal- 
ance, which has a st rong impact on the line force and, consequently, wind acceleration. 



Muijres et al.l (120111) studied the effect of different portions of clumped and unclumped 
parts of the wind and different clump sizes (assuming that clumps are cubes). They 
used the Monte Carlo wind code of Vink et al. (2000, 2001) and found that optically 
thin clumps increase the mass-loss rate, whereas optically thick clumps decrease the 
mass loss rate. 
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5. Conclusions 

The problem of P v mass-loss rate determination has not been satisfactorily solved yet. 
While X-rays do not seem to alter the ionization balance of P v, the case of XUV 
radiation has still to be tested by detailed consistent NLTE calculations. 

The dependence of mass-loss rates on clumping has still not been analysed in a big 
detail yet. Available results show that clumping below the critical point may have strong 
influence on predicted mass-loss rates. Although results of hydrodynamical simulations 
do not support clumping below the sonic point, there is a possibility of creating clumps 
by other mechanisms. Detailed hydrodynamical calculations are needed. 
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